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The majority of magnetic materials possess some degree of magnetic anisotropy, either at the 
level of a single ion, or in the exchange interactions between different magnetic ions. Where these 
exchange interactions are also frustrated, the competition between them and anisotropy can stabilize 
a wide variety of new phases in applied magnetic field. Motivated by the hexagonal delafossite 2H- 
AgNi02, we study the Heisenberg antiferromagnet on a layered triangular lattice with competing 
first- and second-neighbour interactions and single-ion easy-axis anisotropy. Using a combination of 
classical Monte Carlo simulation, mean-field analysis, and Landau theory, we establish the magnetic 
phase diagram of this model as a function of temperature and magnetic field for a fixed ratio of 
exchange interactions, but with values of easy-axis anisotropy D extending from the Heisenberg 
(D—0) to the Ising (D=oo) limits. We uncover a rich variety of different magnetic phases. These 
include several phases which are magnetic supersolids (in the sense of Matsuda and Tstuneto or 
Liu and Fisher), one of which may already have been observed in AgNi02- We explore how this 
particular supersolid arises through the closing of a gap in the spin-wave spectrum, and how it 
competes with rival collinear phases as the easy-axis anisotropy is increased. The finite temperature 
properties of this phase are found to be different from those of any previously studied magnetic 
supersolid. 



PACS numbers: 67.80. kb, 75.10.-b, 75.10.Hk 

I. INTRODUCTION 

The very real possibility of finding new quantum 
phases and excitations has made frustrated magnetism 
one of the most exciting and dynamic research areas in 
contemporary condensed matter physics^l Often, frus- 
trated magnets have the potential to support many differ- 
ent ground states, and can be tuned at will between these 
states by applying a magnetic field. Relative to other 
approaches, such as chemical doping, this field control 
of phase transitions has the great advantage that it does 
not introduce new complications — such as disorder — 
into the system. The advent of powerful new magnets, in 
some cases capable of reaching many hundreds of Tesla, 
has made it possible to study new magnetic phases and 
the (quantum) phase transitions between them in a wide 
range of materials for the first time. These high-field 
experiments have already yielded new insights into phe- 
nomena as diverse as Bose-Einstein condensation^, and 
the delicate interplay between magnetism and structural 
transitions in Cr spinels^. 

Triangular lattice antiferromagnets have been central 
to our understanding of frustration since the pioneer- 
ing work of Wannier^-. A wide range of novel phases 
have first been proposed in the context of the trian- 
gular lattice, including the celebrated "resonating va- 
lence bond" (RVB) spin-liquid ^, tw o-dimensional chiral 
phases^, magnetisation plateauxP^, and a collection of 
different multipolar states^HSl. Recently, spin-liquid be- 
haviour has been established experimentally in a variety 
of triangular-lattice systems, including the quasi-2D or- 
ganics K-(ET) 2 Cu 2 (CN)^and EtMe 3 Sb[Pd(dmit) 2 ] 2 M1 , 
thin films of 3 He [Ref. [T7], and the layered pnic- 



tide NiGa 2 AsJ^. 

While the majority of this work assumes SU{2) invari- 
ance — i.e. the absence of magnetic anisotropy — this 
is rarely true of real materials. Generically, these exhibit 
some degree of magnetic anisotropy, either at the level of 
a single ion, or in their exchange interactions^. It has 
been known for a long time that the interplay between 
magnetic field, frustration and magnetic anisotropy can 
lead to novel magnetic phases on the triangular lattice, 
even in simple classical models. The classical Heisenberg 
antiferromagnet with easy-axis anisot ropy, first stud- 
ied by Miyashita and Kawamur al 20 * 21 !, exhibits a com- 
petition between the coplanarity favoured by magnetic 
interactions^, and collinear order driven by anisotropy 
and magnetic fields. More recent studies have high- 
lighted successive Berezinskii-Kosterlitz -Thou less phase 
transitions as a function of temperatur e 1 23 * 24 ! and exotic 
quantum effects at large single-ion anisotropy^. 

One possible outcome of the competition between ex- 
change interactions, magnetic anisotropy and magnetic 
field is a state which breaks both spin and lattice symme- 
tries, exhibiting a finite (staggered) magnetisation per- 
pendicular to the field, and broken translational sym- 
metry. Matsuda and Tstuneto^ and Liu and Fisher^l 
argued, via a quantum lattice gas mapping, that such 
a mixed-symmetry phase could be an exact analogue of 
the supersolid phase proposed in 4 He [Ref. [251 - f5U] . In the 
model considered by Liu and Fisher, a magnetic super- 
solid was found to interpolate between a collinear anti- 
ferromagnet and its high-field, canted, "spin-flop" state. 
This magnetic supersolid was argued to terminate in a 
finite-temperature tetracritical point, where collinear, su- 
persolid, spin-flop and paramagnetic phases meelP^. 



Recent, controversial, experiments on 4 He [Ref. [3T] 
have lead to a renaissance of interest in supersolid 
phases^ and their magnetic analogues. Quantum Monte 
Carlo simulations now provide clear evidence for mag- 
netic supersolids in S=l/2, easy-axis models on the trian- 
gular lattice 33 '®!. In this case, following the original logic 
of Andreev^l, the transition into a magnetic supersolid 
can be interpreted as a field-driven Bose-Einstein conden- 
sation of magnons within a magnetic solid — a collinear, 
magnetically ordered phase. Magnetic s uperso lids have 
also been associated with biconical phaseipZESl. 

Among triangular lattice magnets, AgNi02 is of special 
interest. Combining metallicity and localised magnetism, 
its ground state is a collinear antiferromagnetic "stripe" 
phased When applying magnetic field a complicated 
and as of yet unexplained cascade of phase transitions 
is observed^. In a previous work we have argued that 
the low-field transition observed in AgNiC>2 converts the 
collinear ground state into a novel state which is a mag- 
netic supersolid in the sense of Matsuda and Tstuneto or 
Liu and Fisher^S. In this case, the expected first-order 
spin-flop transition is replaced by a continuous transi- 
tion into a magnetic supersolid. This occurs through the 
Bose-Einstein condensation of magnons with a finite mo- 
mentum, leading to a novel magnetic supersolid distinct 
from any of the cases considered above. 

In this paper we extend the analysis began in Ref. [41] 
to explore systematically how competing exchange inter- 
actions and easy-axis anisotropy provide a route to novel 
magnetic phases under applied field in, e.g., hexagonal 
delafossites. We consider explicitly the model used to fit 
low-energy magnetic excitations in AgNi02, with com- 
peting first-neighbour (Ji) and second-neighbour (J2) 
exchange interactions on a triangular lattice, small but 
finite inter-layer coupling (Jl), and a varying degree 
of easy-axis anisotropy D [Ref. 02]. While there have 
been some theoretical studies on the interplay of further- 
neighbour interactions and anisotropy on the triangular 
lattice^ES, the possibility of finding new phases in mag- 
netic field remains largely unexplored. 

In this paper we constrain the ratio between exchange 
parameters to values relevant to AgNi0 2 and explore 
how the magnetic phase diagram evolves as a function 
of easy-axis anisotropy. The resulting behaviour is ex- 
tremely rich, exhibiting a wide variety of magnetic super- 
solid and collinear plateau phases, (almost) all of which 
share a common four-site unit cell. Since our main in- 
terest is in finite-temperature phases and phase transi- 
tions, we consider the classical limit S — > 00, which opens 
the problem to large-scale Monte Carlo simulation. This 
is complemented by zero-temperature mean-field theory, 
low-temperature spin-wave expansions, and Landau the- 
ory. Particular attention is paid to the novel magnetic 
supersolid phase introduced in Ref. 1411 This is shown to 
have qualitatively different properties at finite tempera- 
ture from any previously studied supersolid, and to be 
remarkably robust against changes in parameters. 

The plan of the paper is as follows: in Sectionlnlwe in- 
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FIG. 1: (color online) (a) Layered triangular lattice considered in 
this paper, with first- and second- neighbour interactions J\ and 
J2 in each plane, coupled by interactions between perpendicular 
neighbours J±. (b) Second-neighbour interactions transform the 
natural unit cell within the plane into a diamond plaquette. 



troduce the model, discuss the choice of parameters and 
present mean-field arguments. We also briefly outline 
the numerical and analytical methods employed. In Sec- 
tions. |niHlX]we explore the finite-temperature and finite- 
field behaviour of the model for representative values of 
easy-axis anisotropy, interpolating from the Heisenberg 
model [D = 0] to the Ising limit [D = 00]. The main 
results of the paper are summarised in a correspond- 
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D = 0.65], Fig7[28| [D = 1.5] and Fig, _ 
We summarise the key points of the results and discuss 
the remaining open questions in Section [X] In Section [Xl| 
we consider the relevance of these results to experiments 
on AgNi02- We conclude in Section [XII| Technical de- 
tails of calculations of spin stiffness and low-temperature 
spin wave expansions are discussed in Appendix [A] and 
Appendix [S] 



II. MODEL AND METHODS 

We consider the classical Heisenberg model on a lay- 
ered triangular lattice 

n =Ji s i- s i + J ^Y1 s ^ + J ^Y1 s - s j 

{ij)2 (ij)± 

-D^iStf-hY^S?, (1) 

i i 

with competing antiferromagnetic interactions J\ > 
and J2 > on the first and second neighbour bonds of 
a triangular lattice plane [Fig. [I] . Neighbouring planes 
are coupled by an interlayer exchange Jj_. Easy-axis 
anisotropy is characterised by a single-ion term D and 
the magnetic field h is assumed to lie along the z axis. 
We set |S|=S'=1 throughout, and all energy scales in- 
cluding field h and temperature T are measured in units 
of J 1 hereafter. 

We fix first-neighbour interactions to be J\ = 1 and set 
interlayer exchange J±—— 0.15. This leaves three distinct 



free parameters: second- neighbour coupling J 2 , easy- axis 
anisotropy D and external field h. We note that the 
h = 0, J± = 0, D — » 00 limit of Eq. ([!]) corresponds 
to the frustrated Ising model previously studied in e.g. 
Refs.i5]and|lll 

Second-neighbour interactions change the "natural" 
description of the lattice from edge-sharing triangles to 
edge-sharing four-site "diamond" plaquettes, which can 
be thought of as a 2D projection of an irregular tetrahe- 
dron. Thus, at a mean-field level, the exchange terms of 
Eq. can be rewritten in terms of the total magnetisa- 
tion of each plaquette 

T + 7 N/i 

H ex = ^^Y,\ S A + S B + S c + S D \ 2 -4. (2) 



This allows an interesting parallel with the physics of 
tctrahedra embedde d in a pyrochlore lattice in a bilinear- 
biquadratic modeJiJilliSI, explored further below. 

The application of a magnetic field reveals a cascade 
of new states, summarised in Fig. [2j These consist of : 

(a) A two-sublattice, collinear antiferromagnetic "stripe" 
state [Fig. [2ja)], which is the ground state in the 
absence of magnetic field, and which breaks both the 
translational and rotation symmetries of the lattice. 

(b) A two-sublattice canted "spin-flop" state [Fig. [2jd)] 
which breaks both the rotational symmetries of the 
lattice, and spin-rotation symmetry in the S x -S y 
plane. 

(c) A novel magnetic supersolid state [Fig. |2je)] with a 
four-site unit cell which breaks both the translational 
and rotational symmetries of the lattice, and spin- 
rotational symmetry in the S x -S y plane 

(d) A three-sublattice, to=1/3 magnetisation plateau 
[Fig. [2jb)] which breaks the translation symmetries 
of the lattice. 

(e) A four-sublattice, 2:1:1 canted magnetic supersolid 
state [Fig. Etf)], which breaks both the rotational 
symmetries of the lattice and spin-rotational symme- 
try in the S x -S v plane 

(f) A four-sublattice, m=l/2 magnetisation plateau 
[Fig. [2jc)] which breaks both the translational and 
rotational symmetries of the lattice. 

(g) A four-sublattice, 3:1 canted magnetic supersolid 
state [Fig. [2^g)] , which breaks both the rotational 
symmetries of the lattice and spin-rotational symme- 
try in the S x -S v plane 

Each of these phase will be described in more detail be- 
low. In what follows we shall focus in particular on the 
novel magnetic supersolid state [Fig. [2jc)], previously in- 
troduced in Ref. UH as an explanation for the low-field 
transition observed in AgNiOspSl. 
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FIG. 2: (color online) Cartoons of different three- and four- 
sublattice ground states found at mean-field level as a function of 
field, (a) Collinear stripe phase, (b) "Spin-flopped" version of the 
stripe phase, (c) novel magnetic supersolid state, the main focus 
of this paper, (d) m=l/3 magnetisation plateau, (e) 2:1:1 canted 
phase, (f) m=l/2 magnetisation plateau, (g) 3:1 canted phase. 
The staggered moment in the S x —S v plane is shown as a vector 
for states (b), (c), (e) and (g). 

To get a first idea of how these phases fit together, let 
us first consider the phase diagram in the absence of field, 
for a range of J 2 and -D, within zero-temperature mean 
field theory [Fig. [3j. For vanishing magnetic anisotropy 
(D = 0) and J 2 < 1/8, we find the well-known coplanar 
three-sublattice 120° ground state 8 . For 1/8 < J 2 < 1 
this gives way to a set of classically-degenerate two- 
sublattice antiferromagnetic configurations. These are 
locked together, by thermal^ or quant um flu ctuations, 
to form collinear "stripe" order [Figj4^a)j 50 51 . For J 2 >1 
the ground state is a one-dimensional coplanar spiral 
with wave number q = (2/^3) cos" 1 [-(J x + J 2 )/(2J 2 )] 
[Ref. I49l50j . The inclusion of easy-axis anisotropy fur- 
ther stabilizes the collinear stripe state, and converts the 
120° state and coplanar spirals into a three-sublattice 
coplanar farPS and a one-dimensional helicoidal phase, 
respectively [Fig. [3]. 
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FIG. 3: (color online) Zero field J2—D phase diagram for Ji=l 
and J^=-0.15 from zero temperature mean-field theory. The model 
supports three different ground states : a three-sublattice coplanar 
state, a two-sublattice collinear "stripe" antiferromagnet and an 
incommensurate helicoidal phase. 

Now let us consider how this zero-temperature phase 
diagram evolves with field. A reasonable expectation 
would be that, for sufficiently large field, the collinear 
stripe state [Fig. [2ja)] would undergo a first-order 
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FIG. 4: (color online) (a) Collinear stripe phase favoured by 
second-neighbour interactions, showing the three possible stripe 
orientations, (b) Three-sublattice m=l/3 plateau state with two 
spins up and one down, (c) First Brillouin zone in q x -<]y plane 
showing momentum vectors M and M' associated with the stripe 
phases. The blue rectangle shows the first magnetic Brillouin zone 
for the stripe phase with ordering vector M (top cartoon), (d) 
The red hexagon shows the first magnetic Brillouin zone for a 
generic four-sublattice state centred on M'. (e) Ordering vectors 
K of the m=l/3 plateau state and related first magnetic Brillouin 
zone, shaded blue. 



"spin-flop" transition into a two-sublattice canted state 
[Fig. [2jb)] . However this is not what happens. The first 
new phase encountered as magnetic field is increased, 
within mean field theory, is shown as a function of J2 
and D in Fig. [5j At small J2, we find the familiar 
collinear m = 1/3 magnetisation plateau [Fig. [4^c)] and 
for large D, moderate J2, a collinear m — 1/2 magneti- 
sation plateau [Fig. [2jf)]; for large J 2 and small D the 
helicoidal phase interpolates all the way to saturation. 
However, away from the margins of the phase diagram, 
a novel four-sublattice, partially canted phase [Fig. [2^c)] 
completely dominates, displacing the more conventional 
spin-flop phase. This partially canted state is a magnetic 
supersolid phase, according to the definition of Matsuda 
and Tsunetc^or Liu and Fisher^. While the stripe state 
breaks only discrete lattice symmetries, and the canted 
spin-flop state only continuous spin-rotational ones, the 
partially canted supersolid breaks both. The reasons for 
the emergence of this unusual — and very robust — new 
state form the main focus of this work. 

For the purposes of this paper we choose to fix 
the nearest-neighbour interaction J2=0.15 and vary the 
anisotropy strength, selecting values along the vertical 
dashed line in Fig. [3] This line includes the param- 
eters J2=0.15Ji and D=2/3Ji relevant to the magnet 
AgNi0 2 [Ref. Elj discussed abovcP. Moreover, since the 
J2 = 0.15 line lies close to several phase boundaries in 
J-2-D space, it opens the possibility of field-driven com- 
petition between different phases. 

The main results of this paper are a set of finite- 
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FIG. 5: (color online) Finite-field J2—D phase diagram for J\=l 
and Jx=-0.15, showing the first mean-field instability with applied 
field. Three different states are found : m=l/3 and m=l/2 mag- 
netisation plateaux and magnetic supersolid state. Red dots along 
dashed vertical line show studied anisotropy values. 



temperature, finite- field phase diagrams for Eq. ob- 
tained using classical Monte Carlo (MC) simulation. Be- 
cause of the strong magnetic frustration and magnetic 
anisotropy, and a large number of first-order phase tran- 
sitions, MC simulations of this model are very challeng- 
ing. To overcome severe slowing down in simulation dy- 
namics we employ a parallel tempering MC scheme^, 
which combines local Metropolis updates with parallel 
simulation of multiple replicas of the system at differ- 
ent temperatures. The exchange of replicas at different 
temperatures allows the system to avoid local free energy 
minima traps. To further reduce the autocorrelation be- 
tween spin configurations, especially in the equilibration 
phase, every local-update sweep is followed by successive 
over-relaxation sweeps^. These are entirely determinis- 
tic, and comprise the reflection of each spin in turn in 
the plane formed by the S* easy axis and the local field 
from its neighbouring spins. This over-relaxation scheme 
is a micro-canonical update, and reversible, and there- 
fore the global Markov chain for parallel tempering and 
over-relaxation will also obey detailed balance. 

Simulations of from 48 to 160 replicas were performed 
for rhombohedral clusters of 3L x 3L x L = 9L 3 spins 
with periodic boundary conditions, where L=4,6,8,10 is 
the cluster thickness, i.e. counts the number of triangular 
lattice planes. This cluster geometry was chosen to reflect 
the higher correlation length in the x— y plane relative to 
the z axis for Jl/Ji«0.1. Typical simulations involved 
4x 10 6 steps, half of which were discarded for thermaliza- 
tion. Each step consisted of one local-update sweep of the 
lattice followed by two over-relaxation sweeps, with repli- 
cas at different temperatures exchanged every 10 steps. 
Typically, simulations were started from initial random 
configurations. Nevertheless, the competition between 
three and four-sublattice states leads to strongly first or- 
der transitions involving large internal energy disconti- 
nuities at low temperature, difficult to overcome even 
with parallel tempering. In order to overcome this, the 
results in this region were obtained through careful com- 
parison between several runs with different initial config- 
urations. Some runs were initialised with fully ordered 
states and others with mixed phase configurations be- 
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FIG. 6: (color online) Magnetic phases of a layered triangular-lattice antiferromagnet with exchange interactions Ji=l, J2=0.15, 
Jx=— 0.15 and easy-axis anisotropy D=0. (a) Phase diagram as a function of temperature and magnetic field. Open symbols on the 
h-axis show transitions obtained in mean-field theory. Phase boundaries at finite temperature are obtained from Monte Carlo simulation 
for a cluster of 24x24x8 spins, and determined by peaks in the relevant order parameter susceptibilities. All phase transitions are first- 
order, except where shown with a dashed line. All phases share a 4-sublattice unit-cell, (b) Phase diagram as a function of temperature 
and magnetization. Solid lines running left-right show cuts at constant magnetic field h taken from simulations. The coexistence regions 
associated with first order phase transitions are coloured white. 



tween both types of order. 

These simulations are complemented by analytic cal- 
culations at low temperature. The T—0 ground state of 
the model can be calculated within mean-field theory for 
2, 3 and 4 sublattice states. This provides a direct check 
on the T=0 limit of the phase transitions found in sim- 
ulations. The Gaussian fluctuations about these mean- 
field states are calculated using the low-temperature spin- 
wave technique outlined in Appendix A. We use this ap- 
proach to calculate magnetisation curves at finite field, 
and to parametrize a Landau theory for continuous phase 
transitions. These results are compared with phase 
boundaries obtained from simulation in the phase dia- 
grams below. 



III. HEISENBERG LIMIT, D=0 

The two natural limits of the Hamiltonian Eq. ([I]) are 
the Heisenberg (D=Q) and Ising (D=oo) models. We 
consider a range of D interpolating between these, begin- 
ning with a pure Heisenberg model. For D=0, the phase 
diagram is dictated by the order-from-disorder selection 
of a small number of magnetically ordered phases from a 
disordered manifold of competing ground states. In the 



absence of magnetic field, like previous authora^^l, we 
find that fluctuations favour a collinear antiferromagnetic 
"stripe" ground state. Once magnetic field is applied, 
these stripes cant. At a field of ft,=4(J 1 + J 2 ) = 4.6 
(T=0 + ), there is a first order transition from the 2- 
sublattice canted state into a 4-sublattice collinear half- 
magnetization plateau state. As field is increased further, 
this plateau undergoes a continuous phase transition into 
a 3:1 canted state, which interpolates to saturation. Fi- 
nally the system saturates at a field of h=9.2 (T=0). 
These results are summarised in Fig. [6j and described in 
more detail below. 



In the absence of magnetic field (h=0) the zero- 
temperature mean-field ground state of Eq. ( l]) is a man- 
ifold of states built of two decoupled Neel sublattices 
[Fig. Ra)]. Thermal fluctuations select the collinear 
stripe phase [Fig. [7Fb)] which maximizes the entropy ac- 
cessible to spin wave excitations. This stripe state breaks 
both 0(3) spin-rotation symmetry and the (discrete) ro- 
tation symmetry of the lattice — the stripes can be ori- 
entated in three different ways [Fig.ga)]. To study the 
finite temperature transition into this state, we there- 
fore introduce a complex order parameter ip = ipi + iip2 
based on a two-dimensional irreducible representation of 




FIG. 7: (color online) Structure of the ground state of the J1—J2 
Heisenberg model on a triangular lattice for J2/J1 > 1/8. (a) De- 
generate manifold of states for T=0, h=0, built of two independent 
Neel sublattices. (b) Collinear antiferromagnetic stripe phase se- 
lected by thermal fluctuations (h=0). (c) "Spin flop" phase formed 
by direct canting of antiferromagnetic stripes in applied magnetic 
field, (d) Supersolid phase formed from the same two canted sub- 
lattices. 



the C3 = Z3 lattice rotation group 

ipi = ^/q^y E 2Sj.Sj + ^ 1 — Sj.Sj + 5 2 — Si.Si-s 1 s 2> 

V>2 = E Sj.Sj +52 - S i .S i _ Sl _s 2 , 

Here the primitive vectors of the triangular lattice are 
61 = (1,0), 62 — (1/2, v / 3/2). Since parallel tempering 
effectively restores this lattice symmetry, we measure the 
magnitude of the order parameter 



(3) 



which is normalised to 4/\/6 for perfect stripe order. In 
Fig. [H] we present MC simulation results for the finite 
temperature transition from the paramagnetic into the 
collinear antiferromagnetic state. The transition is found 
to be first order, as expected for a 3-state Potts model in 
3D (see e.g. Ref-HU). 

Applying a magnetic field causes the spins in each Neel 
sublattice to cant. This can happen in two distinct ways. 
Both stripes can cant equally, to give the spin flop state 
shown in Fig. [7|c). Alternatively, the sublattices can 
rotate so that half of the spins are aligned with the mag- 
netic field, while the others are canted away from it 
the supersolid state shown in Fig.CTd). These two states 
are degenerate, but thermal fluctuations favour the spin- 
flop state over the supersolid, since it has a higher density 
of spin wave excitations at low energies. 

In the spin-flop phase all spins have equal S z com- 
ponents, and the lattice symmetries are broken only by 
spin components in the S x -S y plane. We therefore mod- 
ify Eq. (|3| to read 



62 ' 



(b) 




0.1 



0.2 
T 



0.3 



0.4 



FIG. 8: (color online) Phase transition at D = 0, h = from the 
paramagnet into the isotropic collinear antiferromagnetic stripe 
state, (a) Stripe order parameter from broken Z3 lattice sym- 
metry, Eq. |3j, and (b) heat capacity. The first-order behaviour is 
clear in the bimodal distribution of the energy histograms close to 
the critical temperature, inset to (a). 



where Sj~ = (Sf, Sf), and consider the order parameter 



Oi_ = 



l^l 



(4) 



The spin-flop phase also separately breaks spin-rotation 
symmetry in the S x -S v plane. This can be determined 
using an order parameter for the staggered in-plane mag- 
netization 



O 



U(l) 



1 

N 



N/4 

E 





s- 1 



si I 



(5) 



where A and B label the two sublattices of the spin-flop 
state, as illustrated in Fig. [2] (b) . The breaking of spin- 
rotation symmetry implies the existence of a finite spin 
stiffness p$, as defined in Appendix [A} 

At higher values of magnetic field, thermal fluctua- 
tions select another new ordered state from a large set of 
competing alternatives. This is a collinear state with a 
4-site unit cell in which three spins point "up" and one 
"down" — commonly referred to as a half-magnetization 
(m=l/2) plateau [Fig. [2](f )] . In the limit of vanish- 
ing temperature, this first occurs at a magnetic field 
value of h=4(Ji + J2)=4.6. This collinear phase breaks 
translational symmetries, with associated ordering vec- 
tors {qMj-, cf. Fig. [3J The S z -S z component of the 
static spin structure factor at these wave vectors 



S zz (q M ) = 



E 

{qm} 



1 

N 



E^ 



(6) 
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FIG. 9: (color online) Selection of simulation results for D = 0, 
h = 4.0, spanning the canted "spin-flop" phase, collinear m=l/2 
plateau, 3:1 canted state, and high temperature paramagnet. (a) 
Z3 order parameter Oz 3 associated with broken lattice rotation 
symmetry in the spin flop state [Eq. pi], and spin structure factor 
cS 2z (qM) associated with broken translational symmetry in the 3:1 
canted and collinear m=l/2 plateau states [Eq. ijfjj]; (b) in-plane 
staggered magnetization [Eq. 151]; (c) spin stiffness pg; (d) 

heat capacity C^; (e) magnetic susceptibility \- The inset to (a) 
shows the energy distribution at the transition from paramagnet 
to 3:1 canted phase. The inset to (b) shows the crossing of the 
Binder cumulants associated with Oy(i) at the transition from the 
3:1 canted phase to the collinear m=l/2 plateau state. 



acts as (the square of) an order parameter for this state. 
Since the state is collinear, and aligned with the magnetic 



field, its spin stiffness ps = 0. We note that, since the 
"down" spin is unique within the 4-site unit cell, this 
collinear state also breaks the permutation symmetries 
of the bonds within the unit cell. It is therefore possible 
to define a complementary Z4 order parameter for this 
state, based on the irreps of the relevant permutation 
group. For present purposes, however, the combination 
of a vanishing ps and a finite S zz (qfj\) are sufficient to 
identify the state. 

The catalogue of ordered phases for 23=0 is completed 
by a 3:1 canted phase [Fig. [2jg)], which is connected to 
collinear m=l/2 plateau by a continuous phase transi- 
tion as a function of magnetic field. All spin wave modes 
are gapped in the collinear plateau state. With increas- 
ing magnetic field, the gap associated with one of the 
zone-centre spin waves closes. Since this spin wave is a 
transverse excitation, it converts the collinear state into 
one with a staggered magnetization in the S x -S y plane, 
and a finite spin stiffness. And, since, it is a zone-centre 
mode, the 3:1 canted phase inherits the broken transla- 
tional symmetry of the m=l/2 plateau. The resulting 
state is therefore a supersolid, in the sense of Matsuda 
and Tstuneto or Liu and Fisher, breaking both transla- 
tional and spin-rotational symmetries. 

In Fig. [9] we show a selection of finite temperature sim- 
ulation results for h—4.0, a cut through the T-h plane 
which spans all three ordered states. The transition from 
the high temperature paramagnet to the 3:1 canted phase 
at T=0. 290(3) can be seen in the temperature depen- 
dence of S zz (c[ M ) [Fig. (9ta)], staggered magnetization 
[Fig.[9](b)], spin stiffness ps [Fig. [9^c)] , and a sharp 
feature in <Jh [Fig- [^d)]. While symmetry permits a 
second-order phase transition, the phase transition ob- 
served in simulation for this value of h is first-order, as 
shown by the energy histograms in the inset to Fig. [9ja) . 
The transition from 3:1 canted phase to m=l/2 plateau is 
distinguished by a collapse of the spin stiffness [Fig.[9jc)], 
0(7(1) [Fig- [Sp 3 )] ancl a suppression of the magnetic sus- 
ceptibility x [Fig- (2]l e )]- The transition is continuous, 
with the Binder cumulant for the two-component U(X) 

order parameter t/4 = 1 — „ lr !i (1 \ 2 crossing at a single 

temperature T=0. 210(3) for a range of system sizes [in- 
set to Fig. [9jb)]. A strong finite size dependence of this 
Binder cumulant is observed for all observed continuous 
phase transitions. These finite-size corrections are espe- 
cially large in ordered phases (e.g. plateau) where the 
relevant U(l) order parameter vanishes in the thermody- 
namic limit. 

Finally, the abrupt, first-order transition from the 
m=l/2 plateau to the spin flop phase at T—0. 100(3) 
is visible in the sudden onset of 0% 3 [Fig. |9ja)], 0£/(i) 
[Fig.KMb)] and spin stiffness [Fig.[9^c)]. Here a continuous 
phase transition can be ruled out on symmetry grounds. 

The magnetic phase diagram can also be plotted as a 
function of temperature and magnetisation [Fig.^b)]. In 
this case a discontinuity in the magnetisation as a func- 
tion of temperature or magnetic field implies a first-order 
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FIG. 10: (color online) Magnetic phases of a layered triangular-lattice antiferromagnet with exchange interactions Ji=l, J2=0.15, 
Jx=— 0.15 and easy-axis anisotropy D=0.02. (a) Phase diagram as a function of temperature and magnetic field. Open symbols on the 
h-axis show transitions obtained in mean-field theory. Phase boundaries at finite temperature are obtained from Monte Carlo simulation for 
a cluster of 24x24x8 spins, and determined by peaks in the relevant order parameter susceptibilities. All phase transitions are first-order, 
except where shown with a dashed line. Thick purple dashed line is obtained through a Landau expansion for the supersolid transition, 
(b) Phase diagram as a function of temperature and magnetization. Solid lines running left-right show cuts at constant magnetic field h 
taken from simulations. Thick purple dashed lines show phase boundaries obtained through a low-T expansion. Finite anisotropy lifts 
T=0 degeneracy in favour of the stripe and supersolid phases, relegating the spin-flop to finite-temperature only. 



phase transition, and leads to regions of phase coexis- 
tence in the T-m plane, coloured white in Fig. pi b). The 
transition from the paramagnet to the m=l/2 plateau 
is thus seen to be first order, although with a very small 
magnetisation jump. The phase transition from the spin- 
flop phase to the paramagnet is also first order for h>0. 

The transition between the m=l/2 plateau and the 3:1 
canted state is continuous, as is the finite-temperature 
transition from the 3:1 canted state to the paramagnet, 
for h >5 (m >0.5). The 3:1 canted phase finally saturates 
at zero temperature, at a field h$AT — 9.2. 

It is interesting to note that the same succession of 
phases in magnetic fiel d occ urs in a model of spin-lattice 
coupling in Cr spinels^- 47 * 48 !. In that case the four-site cell 
is the tetrahedron from the which the pyrochlore lattice 
is built, and the role of thermal fluctuations is played by 
a biquadratic spin interaction. 



IV. SMALL ANISOTROPY, D=0.02 

The phase diagram for D—0.02, presented in Fig. [Toj 
is extremely rich. It exhibits all but one of the phases 
catalogued in Fig. [2] — a collinear stripe phase at low 



values of field, giving way to either a supersolid or a 
spin-flop phase as magnetic field is increased; a novel 
2:1:1 canted state, a collinear m=l/2 plateau state and 
finally, approaching saturation, a 3:1 canted phase. The 
key to understanding the richness of this magnetic phase 
diagram is to recognise that none of these phases found 
for D = are selected by energetic considerations alone 
— all of them owe their stability to thermal fluctuations. 
The difference in entropy between these phases and the 
other degenerate alternatives is, however, very small. For 
this reason the introduction of even a vanishingly small 
easy-axis anisotropy has profound consequences for the 
phase diagram, especially at low temperatures. 

To explore how this works, we first consider the T=0, 
mean-field energies per spin of the two competing states 
at low magnetic field, the canted spin flop state (SPF) 
and the partially canted magnetic supersolid (SSD) 

Zs PF = -Ji-Ji + J,- m{Ji+ h j 2) „ AD , (7) 

= -1 6 ( J i + J 2 ) ( J i + J 2 - J± ) - 8 D ( Ji + J 2 + J± ) + 4 D 2 
SSD 16( Ji + J 2 ) - 8D 



16(Ji + J 2 )-8D' 
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Both of these state are built of the same, canted Neel 
sublattices, but with different relative orientations, as il- 
lustrated in Fig. [7] (c) and (d). Their energies should 
be compared with the (field-independent) energy of the 
collincar stripe state (STR) 



E S 



TR 



- Ji - J 2 + J± - D. 



(9) 



In the absence of anisotropy (i.e. for D=0), the spin 
flop and supersolid states are degenerate. It is thermal 
fluctuations which select both the collinear stripe state in 
zero field, and the canted spin flop state for h>0. Intro- 
ducing a finite anisotropy D singles out the stripe phase 
as the ground state, and lifts the degeneracy between the 
canted spin flop and supersolid states. However, while 
thermal fluctuations favour the spin flop state, anisotropy 
favours the supersolid, since the spins lie closer to the 
easy axis. This sets up an interesting tension between 
energy and entropy which, for very small values of D, 
compete on an equal footing. 



Q 
lD 




FIG. 11: (color online) Comparison of zero-tempcrature mean 
field energies for the canted spin flop (blue, dashed line) and mag- 
netic supersolid (red solid line) states for parameters J2=0.15 and 
Jx=-0.15 as a function of magnetic field h. A continuous phase 
transition into the supersolid state occurs at hsSD=2D, anticipat- 
ing the first order spin flop transition which would otherwise have 
occurred at /ispf, Eq. ( |11[ |. 



degeneracy is lifted by an applied field in favour of M' , 
with the spin gaps evolving as 

A M = D + 2(Ji + Ji) - i>/l6(Ji + J2) 2 + h 2 , (12) 

A M ,=D-^. (13) 

The gap at M' closes at a field hsso=2D, while the gap 
at M would, hypothetically, close at /ispf [Eq. (11)]. It 
is therefore the closing of the gap at M' which medi- 
ates the phase transition. And, since this spin wave is a 
transverse mode with finite momentum, it converts the 
collinear stripe phase into a supersolid state with a fi- 
nite (staggered) magnetization in the S x -S y plane. The 
spin- wave dispersion at h = ft-sSDi when the gap first 
closes, is illustrated in Fig. [l2j An exactly analogous 
phase transition occurs in the quantum model studied 
with (linear) spin-wave analysis, where it is interpreted 
as Bose-Einstein condensation of magnons 41 , to give a 
condensate which breaks the transverse U(l) symmetry. 
Similar mode-softening transitions have been observed in 
other models of frustrated triangular lattices^. 




FIG. 12: (color online) Classical spin- wave dispersion (S x modes 
only) for J2=0.15, Jj_=-0.15 and D = 0.02, at a magnetic field 
h=2D associated with continuous transition from collinear stripe 
phase into supersolid through closing of the spin gap at M'. Inset 
shows low-energy detail. 



At zero temperature phase transitions are controlled 
by energy alone. The first transition as a function of 
magnetic field is into the supersolid state, at a critical 
field of 

^ssd = 2D, (10) 

In contrast, the field at which the spin flop state first 
becomes energetically favourable is 

h SPF = 4\/ D(Ji + J 2 - D/A). (11) 

Moreover, the transition into the supersolid is continu- 
ous, whereas the (avoided) spin flop transition is first or- 
der. This hierarchy of transitions is illustrated in Fig.[TT| 
The seeds of this supersolid state can also be found 
in a spin-wave expansion around the collinear state for 
D>0. At zero field this yields a gapped dispersion with 
minima at both the ordering momentum vector M and 
at the symmetry related momentum M' [FigQb)]. This 



At finite temperature entropy again enters into the ar- 
gument. Entropy favours both the collinear stripe state 
and the canted spin flop state over the supersolid. This 
has two consequences. At low temperatures, it leads to a 
critical field /issd(T') which increases with temperature. 
And, at higher temperatures, for sufficiently small D, 
entropy can drive a first order transition from the super- 
solid into the spin flop state. We can describe the first of 
these two effects quantitatively within a Landau theory. 

In analogy with Eq. §5§, we define an order parameter 

N/4 


where C and D label the two canted sublattices of the 
supersolid state, as illustrated in Fig. ^c). The Landau 
expansion of the free energy can then be written as 

J r = J r o + ~\O u(l) \ 2 + \\O u{1) \ i + ... , (15) 
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FIG. 13: (color online) Magnetic phases of a layered triangular-lattice antiferromagnet with exchange interactions Ji=l, J2=0.15, 
Jx=— 0.15 and easy-axis anisotropy D=0.25. (a) Phase diagram as a function of temperature and magnetic field. Open symbols on the 
h-axis show transitions obtained in mean-field theory. Phase boundaries at finite temperature are obtained from Monte Carlo simulation for 
a cluster of 24x24x8 spins, and determined by peaks in the relevant order parameter susceptibilities. All phase transitions are first-order, 
except where shown with a dashed line. Thick purple dashed line is obtained through a Landau expansion for the supersolid transition, 
(b) Phase diagram as a function of temperature and magnetization. Solid lines running left-right show cuts at constant magnetic field h 
taken from simulations. The coexistence regions associated with first order phase transitions are coloured white. Thick purple dashed lines 
show phase boundaries obtained through a low-T expansion. Easy-axis anisotropy stabilizes the new phases for a wide-range of values. 



where the coefficient a — hsso(T) — h is calculated exactly 
within a spin- wave expansion about the supersolid state. 
Details of this are given in Appendix [B] For D=0.02 we 
find a critical field of 



hsso(T) = 0.04 + 3.02(4) x T. 



(16) 



The increase of the critical field with temperature is a 
direct consequence of the fact that the supersolid state 
has a lower entropy than the collinear stripe phase. The 
decrease in entropy on entering the supersolid state can 
be traced back to a shift of spectral weight to higher 
energies in the spin-wave spectrum. 

We can also independently estimate hssoiT) from the 
peak in the related order parameter susceptibility 



Xo vm = N 



(0 2 u{1) ) - (Ouw? 



T 



(17) 



which we calculate from Monte Carlo simulation. For 
T<0.2 the analytic and numerical calculations are in 
perfect agreement, as shown by the purple dashed line 
in Fig. 10 a). However for T>0.2 the continuous phase 
transition into the supersolid phase is replaced by a first- 
order transition into the entropically favoured, spin flop 
state. 



As a direct result, there is now a finite temperature 
transition from the supersolid into the spin flop state for 
T«0.2. As might be expected, this phase transition is 
strongly first order, and exhibits marked hysteresis. The 
character of this phase transition is most clearly seen in 
the coexistence region for Tw0.2, mwO.l in Fig. 10 'b). 



In conclusion, although the spin flop state is entirely 
eliminated at low temperatures by the introduction of 
anisotropy, its greater entropy enables it to survive at 
finite temperatures — for sufficiently small values of D. 
From simulations, we estimate that the spin flop phase 
survives at finite temperature for D < 0.045. The great 
wealth of other phases shown in Fig. [10] are also present 
at higher values of D, where simulations are easier to per- 
form. We therefore defer the detailed analysis of these 
phases to the remaining sections of this paper. 



V. MODERATE ANISOTROPY, D=0.25 

The D=0 results showcase the canonical features of 
frustrated magnets in applied magnetic field — see e.g. 
Refs.miandHUl — a collinear plateau state is stabilized at 
intermediate fields by thermal fluctuations, while copla- 
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FIG. 14: (color online) Double transition from the paramagnet 
into the stripe collinear antiferromagnet and then into the super- 
solid phase for D = 0.25, h = 1.0. Translational iS zz (qivi) and 
rotational Z3 [Eq. |18] symmetries are broken at the same tempera- 
ture TssO.37 (a), where system transitions into the collinear phase. 
The first-order character of the outer transition is confirmed by 
the double peak in energy histograms near the critical temperature 
in the inset to (a). A further continuous symmetry is broken at 
T ss 0.16 at the onset of the supersolid phase, as observed in the 
(7(1) order parameter (b) and spin stiffness (c). The different na- 
ture of these transitions is clear in heat capacity (d) and magnetic 
susceptibility (e). 



nar canted phases are found at both lower and higher 
fields. However, in this model even a very small easy- 
anisotropy leads to strikingly different results, as demon- 



strated for Z?=0.02. We now analyse thoroughly the 
phases driven by anisotropy for a representative, mod- 
erately small value of anisotropy £>=0.25. Our results 
are summarized in the phase diagrams shown in Fig. |13| 
For this value of D, the spin flop phase is entirely sup- 
pressed, and the low temperature physics is dictated by 
considerations of energy, rather than entropy. As a func- 
tion of increasing magnetic field, we find a collinear stripe 
state, the supersolid phase described above, an exotic 
2:1:1 canted state, a collinear m=l/2 plateau and finally 
a 3:1 canted state which interpolates to saturation. This 
large set of phases gives rise to a correspondingly large 
number of phase transitions, which we analyze below. 

We begin by considering the finite-temperature phase 
transition from the paramagnet into the collinear stripe 
state at low values of magnetic field, and the subsequent 
low temperature transition from the stripe phase into the 
supersolid. Both of these transition may have been ob- 
served in AgNiC>2 [Ref.|3H]- In Fig. 14 we show a selection 
of simulation results for h — 1.0, and a range of tem- 
peratures spanning all three phases. The transition from 
paramagnet to collinear stripe phase at T=0. 379(4) is ev- 
ident in sharp peaks in both the heat capacity [Fig.|l4^d)] 
and the magnetic susceptibility [Fig. [l4^e)] . This transi- 
tion is first order, as evidenced by a bimodal energy his- 
togram [inset to Fig. 14Ta)]. The collinear stripe phase 
breaks both the translational and the rotational symme- 
try of the lattice. Both of these symmetries are broken 
at the same temperature, as demonstrated in Fig. |14| ^a) 
by the temperature dependence of the spin structure fac- 
tor S zz (qy\) and a suitably modified order parameter for 
lattice rotations 



m 2 



(18) 



where 



. Z ^ \ S rt QZ QZ _ QZ QZ _ QZ QZ 

~ y/6N - D i+i5l ~ D i - D i+S2 ~ D i ■ D i-S 1 -8 2 1 



\ ^ QZ QZ 

/2 N 2^ a i 



Si ■Si-8 1 -8 2 - 



The first order behaviour of the phase transition seen 
here chimes with the known first order phase transition 
from the paramagnet into the collinear stripe phase of 
AgNi0 2 [Ref.HO]. 

The second, internal, phase transition at T=0. 171(4) 
from collinear stripe into the supersolid state is contin- 
uous, exhibiting only shows weak anomalies in specific 
heat [Fig. [l4|d)] and magnetic susceptibility [Fig. 14 [e)]. 
The new broken symmetry — a staggered in-plane mag- 
netization — is heralded by the smooth rise of both the 
(7(1) order parameter 0£/(i) [Fig. 14 ^b)] and the spin 
stiffness [Fig. [l4|c)]. In this context, the susceptibility 



associated with the U(l) order parameter can be anal- 
ysed with a finite-size scaling ansatz 



XOaw - L^xiL^t), 



(19) 
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FIG. 15: (color online) Continuous behaviour of the collinear- 
supersolid transition for D = 0.25, h = 1.0 at Tm0.17. (a) Data 
collapse for the supersolid U(l) order parameter susceptibility us- 
ing 3D XY universality class critical exponents, (b) The Binder 
cumulant for the same order parameter calculated for different lat- 
tice sizes all cross at a single temperature T =0.168(1). 



where t is the reduced temperature t = (T — T c )/T c . 
The critical exponents associated with the order param- 
eter susceptibility (7), and correlation length (v) are ex- 
tracted by fitting this expression to simulation results 
for different lattice sizes [Fig. [15] . The exponents ob- 
tained, ^=1.32(2) and 7=0.67(2}7are in excellent agree- 
ment with the 3D XY universality class. 




FIG. 16: (color online) Classical spin-wave dispersion (S x modes 
only) at _D=0.25 for values of magnetic field associated with contin- 
uous transitions at high field, (a) Closing of the spin gap at M for 
h=1.95, associated with transition from the supersolid phase into 
the 2:1:1 phase, (b) Closing of the spin gap at T for h=5.48, asso- 
ciated with transition from the m=l/2 plateau into the 3:1 canted 
phase. The points M, M' and T are defined in Fig. |ijb). 

In previously studied models with XXZ anisotropy, su- 
persolid phases were found as intermediate phases in- 
terpolatin g betw een Neel antiferromagnetic and spin- 
flop phases^^. However, in the present model, a dif- 
ferent phase emerges above the supersolid. This is a 
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FIG. 17: (color online) Double phase transition at 
D = 0.25, h = 2.5, from the paramagnet into m=l/2 plateau 
and then into 2:1:1 canted state at lower temperature. (a) 
Translational symmetry measured by cS zz (qivi) is broken at TRi0.4 
in the first-order transition - cf. bimodal energy distribution 
(inset) - into the plateau, (b) U(l) continuous symmetry and (c) 
spin stiffness are broken in the continuous transition - cf. crossing 
of U(l) order parameter Binder cumulants in inset to (b) - from 
the plateau into the canted phase at Tft;0.21. (d) The 2:1:1 canted 
state is also characterised by in-plane magnetisation m± (d). The 
different nature of these transitions is resolved in heat capacity (e) 
and magnetic susceptibility (f). 
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2:1:1 canted phase with a four-spin unit cell where two 
spins are parallel to each other and have positive S z , 
while a third is orientated in the negative S z direc- 
tion, and the remaining spin rotates between those po- 
sitions [Fig. |2^e)]. At T=0 this transition is continuous, 
and is observed in both mean-field and spin-wave cal- 
culations, which reveal a soft mode at momentum M 
[Fig. [l6[a)]. With further increase in magnetic field, the 
2:1:1 canted phase evolves smoothly into to the collinear 
m=l/2 plateau. 

Since the 2:1:1 canted phase breaks both translational 
and in-plane symmetries it can also be labelled as a super- 
solid. Furthermore, it possesses a (very small) moment 



(20) 



in the S x —S v plane. This behaviour distinguishes it from 
the lower-field supersolid state, but is also seen in the 
easy-axis nearest-neighbour triangular lattice, where it 
was attributed to the non-trivial degeneracy of the T—0 
ground stated In addition to m±, this phase is charac- 
terised by finite values of the staggered in plane magne- 
tization, spin stiffness ps and structure factor <S Z2 (qM)- 



In Fig. 17 we present simulation results for h — 2.5, 
for a range of temperatures spanning the paramagnetic, 
m=l/2 plateau and 2:1:1 canted phases. The phase tran- 
sition from the paramagnet to the rn=l/2 plateau is 
marked by an abrupt rise in S zz (qM) at T= 0.397(4). 
The transition is first order, as evidenced by a bimodal 
energy histogram [inset to Fig. [TTJa)], and accompa- 
nied by a sharp feature in the heat capacity [Fig. [l7|e)]. 
The phase transition from the m=l/2 plateau to the 
2:1:1 canted state is continuous, with the Binder cumu- 
lants for the associated U(l) order parameter crossing at 
T=0. 211(1) in inset to Fig.^b). 

This continuous phase transition is mediated by a soft 
spin wave mode within the collinear m=l/2 plateau, just 
as the transition into the supersolid is mediated by a 
soft spin-wave mode within the collinear stripe phase. In 
this case the relevant spin-wave gap occurs at the zone 
centre, and closes with decreasing magnetic field, at a 
T=0 critical field of 



h c = 4(J X + J 2 ) - 2D. 



(21) 



Easy-axis anisotropy stabilizes the plateau at T=0 for a 
finite range of field, and its collinearity ensures it is en- 
tropically favoured against the 2:1:1 canted state at finite 
temperature. The critical field therefore slopes down- 
ward with increasing temperature — the opposite of what 
is seen in the transition from the collinear stripe phase 



into the supersolid [Fig. 13 



The transition from the collinear m=l/2 plateau into 
the 3:1 canted state at high field also occurs through 
the condensation of a zone-centre spin wave mode [Fig. 
16 b)], this time at a T=0 critical field of 



h c = 2(Ji + J 2 ) + 2 V / J D 2 +4D(J 1 + J 2 



(Ji + J 2 ) 2 - 



(22) 
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FIG. 18: (color online) Double phase transition at 
D = 0.25, h = 5.3 from the paramagnet to 3:1 canted state 
and then into m=l/2 plateau. The 3:1 canted state is heralded 
by its discrete S zz ((\u\) (a) and continuous U(l) order parameters 
(b), and also by finite spin stiffness (c). This state also has 
a finite in-plane magnetisation (d). The upper transition 

first-order character is clear in the double-peaked energy energy 
distribution close to the critical temperature in inset to (a), while 
the continuous character of the lower transition is shown by the 
Binder cumulant for the U(l) order parameter, in inset to (b). The 
different nature of these transitions is resolved in heat capacity (e) 
and magnetic susceptibility (f). 
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is indeed continuous [Fig. 19 . At T=0 it is easy to see 
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FIG. 19: (color online) Temperature dependence of order param- 
eters associated with 3:1 canted phase close to saturation. The 
continuous nature of the phase transition is confirmed by the single- 
valued energy histogram close to the critical temperature in inset 
to (a), and crossing of Binder cumulants associated with 17(1) order 
parameter in inset to (b). 



This critical field is only weakly dependent on tem- 
perature and, perhaps surprisingly, slopes downwards 
[Fig. ~ 
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As noted in Section III the 3:1 canted phase 



breaks both discrete translational and continuous spin 



rotational symmetries [cf. Fig. 18 . It is therefore a third 
supersolid, in the sense of Matsuda and Tstuneto or Liu 
and Fisher. However, in contrast to the D=0 case, for 
finite anisotropy the 3:1 canted phase possesses a finite 
value of in-plane magnetisation m j_ . 

The 3:1 canted state is the only phase with super- 
solid character which is directly connected to the para- 
magnetic region. This phase transition is clearly first- 
order at low field, e.g. h—5.3, from simulation results 
[Fig. 18 . Nevertheless a continuous phase transition is 



permitted by symmetry, and the transition at high fields 




FIG. 20: (color online) Vanishing discontinuities in magnetisation 
for D=0.25, extracted from temperature cuts in Fig. |13| indicate 
a crossover from first to second-order phase transition at the para- 
magnet - 3:1 canted state transition, where the critical cndpoint is 
found at TsiO.14, h Ri6.5. Line is a guide to the eye. 



that this continuous transition corresponds to the open- 
ing of a spin-wave gap at the 4-sublattice ordering vec- 
tor(s) {qivi}, within the saturated state. For T=0, this 
occurs at 



h S AT = 8(Ji + J 2 ) - 2D. 



(23) 



and the critical field slopes sharply downwards, as re- 
quired by the higher entropy of the paramagnetic phase. 
We conclude that fluctuations drive this continuous phase 
transition first order for T>0.12 (h<6.5) — the point at 
which a finite jump in the magnetization Am z is first 
observed when going from the 3:1 canted phase to para- 



magnet [Fig. 20 



Up to this point we have chiefly concentrated on sym- 
metry breaking at low temperatures, where simulation 
results can be reliably compared with mean-field the- 
ory and spin-wave calculations. We now turn to the 
more delicate question of how the full symmetry of the 
system is recovered with increasing temperature, start- 
ing from the low-field supersolid phase. In the original 
finite-temperature scenario of Liu and Fisher, the su- 
persolid phase extends all the way to the paramagnet, 
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FIG. 21: (color online) Region where four phases compete at in- 
termediate field for D=0.25 (a) and D=0.325 (b) for system size 
L=8. Phase transitions are labelled according to the symmetries 
broken. Three distinct critical points are identified in a narrow field 
range, however these never merge into a single tetracritical point. 
This structure is insensitive to changes in the easy-axis anisotropy 
as long as the T=0 mean field analysis yields the same phases. All 
phase transitions are first-order, except where shown with a dashed 
line. Dashed horizontal lines in (a) show cuts at fixed field studied 
in Fig. [22] 
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FIG. 22: (color online) Successive phase transitions at D=0.25 
for system size L=8 as a function of field as resolved in the sus- 
ceptibilities associated with breakdown of translational symmetry 
cS 2z (qM), spin-rotation symmetry U(l), lattice rotation symmetry 
and in-plane magnetisation mj_. (a) At high field, h=1.85, the 
m=l/2 plateau transforms into the 2:1:1 canted state, followed by 
the supersolid phase, (b) For intermediate field, ?t=1.70, there is a 
direct first-order phase transition from the m=l/2 plateau into the 
supersolid state, as evidenced by the coincidence the susceptibility 
peaks at T rj 0.28. (c) At lower field, /i=1.55, the system transi- 
tions from paramagnct into m=l/2 plateau, then into the stripe 
phase and finally into into supersolid state. Phase transitions are 
first order except where indicated with dashed line. 



where it terminates in a tetracritical poinlP^. Exactly 
at this tetracritical point four phases (collinear antiferro- 
magnet, supersolid, spin-flop and paramagnet) meet, and 
all phase transitions between them are continuous. The 
behaviour of the present model at high temperatures is 
markedly different. From the phase diagrams shown in 
Fig. [l3j it is clear that the collinear phases completely 
engulf the non-collinear ones at high temperature due 
to their higher entropy, so the supersolid phase never 
touches the paramagnet. However, for intermediate val- 
ues of field four phases — the collinear stripe phase, su- 
persolid, 2:1:1 canted phase, and the collinear m—1/2 
plateau — are found in very close proximity to one an- 
other. It is therefore worth asking whether any other 
multicritical points arise in this model. 

In Fig. 21 'a), we present a detailed study of the in- 



termediate field region, 1.2</i<2. Associated field cuts 
for D=0.25 are shown in Fig. [22] The Z3 transition be- 
tween the two collinear states is clearly first order, while 
the transition between the supersolid and 2:1:1 canted 
phases is continuous. The phase transition between the 
collinear stripe state and the supersolid terminates on 
the Z 3 line. The transition between the 2:1:1 canted 
state and the collinear m=l/2 plateau is continuous at 
high fields, but becomes first order shortly before termi- 
nating on the Z 3 line for ft=1.76(l) [Fig. [23]. For /i=1.70 



[Fig. |22[b)] , there is a clear first order transition between 
the m=l/2 plateau and supersolid, where all observed 
susceptibility peaks merge. Given the difficulty in sim- 
ulating this parameter region with many neighbouring 
phases, the possibility of a vanishingly narrow "strip" of 
2:1:1 canted phase extending between the supersolid and 
plateau phases is hard to rule out definitely. However, we 
find no evidence of finite m±, characteristic of the 2:1:1 
state, in this range of temperature and field. 

From this analysis we conclude that the majority of 
these phase transitions remain first order where the four 
phases converge, and that no more than two phases meet 
at a single point via a continuous phase transition. Even 
if the 2:1:1 canted phase were to stretch to lower fields, 
the first-order nature of the phase transitions between 
the m=l/2 plateau and 2:1:1 canted phases, and m=l/2 
plateau and collinear stripe phases would preclude a 
tetracritical point — at most we would have a "bicritical 
endpoint", where two continuous phase transitions meet 
two first-order ones. This behaviour persists throughout 
the anisotropy range where these four phases are present, 
as illustrated in Fig.^b) for D=0.325. 

In the same spirit, it is worth re-examining the fi- 
nite temperature transition between the paramagnet and 
collinear m=l/2 plateau. While generically first-order, 
there exist two points in the phase diagram, /iwl.25 and 
/i«3.2, where the jump in magnetisation between the two 
phases vanishes, indicating the possibility of a continu- 




FIG. 23: (color online) Discontinuities in magnetisation Am z on 
entering the m=l/2 plateau from the supersolid and 2:1:1 canted 
phases for D=0.25 and intermediate field, extracted from temper- 
ature cuts in Fig. |13| The collapse in Ara z for h=1.81(2) indi- 
cates that the transition from the 2:1:1 canted phase to the m=l/2 
plateau becomes first order before merging with the transition from 
the supersolid to the m=l/2 plateau. This precludes a multicritical 
point. The line is a guide to the eye. 
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FIG. 24: (color online) Magnetic phases of a layered triangular-lattice antiferromagnet with Ji=l, J2=0.15, J±=— 0.15 and 
easy-axis anisotropy D=0.5. (a) Phase diagram as a function of temperature and magnetic field. Open symbols on the h-axis show 
transitions obtained in mean-field theory. Phase boundaries at finite temperature are obtained from Monte Carlo simulation for a cluster 
of 24x24x8 spins, and determined by peaks in the relevant order parameter susceptibilities. All phase transitions are first-order, except 
where shown with a dashed line. Thick purple dashed line is obtained through a Landau expansion for the supersolid transition, (b) Phase 
diagram as a function of temperature and magnetization. Solid lines running left-right show cuts at constant magnetic field h taken from 
simulations. The coexistence regions associated with first order phase transitions are coloured white. Thick purple dashed lines show 
phase boundaries obtained through a low-T expansion. Increasing anisotropy gives rise to an m=l/3 plateau with a different unit-cell 
from the other phases. 



ous behaviour [Fig.[l3[b)]. These two points are a robust 
feature of the m=lj2 plateau for all values of D inves- 
tigated. One possible scenario for this behaviour would 
be the emergence of a higher symmetry at these special 
points. Sadly, however, this attractive scenario does not 
survive a closer examination of the data. 

The first point at which Am z —0, at /iwl.25, appears 
to be a critical end point, where the first order transi- 
tion between the collinear m=l/2 plateau and the para- 
magnet terminates on the first order transition into the 
collinear stripe phase. This critical end point is distin- 
guished by a unimodal energy histogram. The second 
point at which Am z —0 occurs for /i«3.20, and can be 
understood in terms of the two different spin wave exci- 
tations which connect the collinear m—1/2 plateau with 
the 2:1:1 canted at lower field, and with the 3:1 canted 
state at higher field. The spin wave-excitation associated 
with the 2:1:1 canted state lowers the magnetization of 
the system, while the excitation associated with the 3:1 
canted state raises the magnetization of the system. 

These spin waves will also determine the sign of the 
magnetization jump at the finite temperature transition 
from the collinear m=l/2 plateau into the paramagnet. 
It follows that there will exist a value of magnetic field 



at which this jump Am 2 changes sign, without any dy- 
namically generated symmetry entering into the problem. 
A close examination of energy histograms and the spin 
structure factor <S(q) at the transition bears out this in- 
terpretation. 

VI. INTERMEDIATE ANISOTROPY, D=0.5 

In order to investigate the robustness of the D—0.25 
results we increase the anisotropy, focusing on a repre- 
sentative value D=0.5. As a consequence a new collinear 
phase emerges at intermediate fields. This is the famil- 
iar three-sublattice state with two spins "up" and one 
"down" [Fig. |2jd)], i.e. a collinear one-third magneti- 
sation plateau with ordering vectors at the corners of 
the Brillouin Zone [Fig[4jd)]. In this case the system 
sacrifices exchange energy associated with the second- 
neighbour interaction to increase its collinearity. 

This m=l/3 plateau first appears within the 2:1:1 
canted phase for Z)>0.33, splitting the canted phase 
in two. Increasing anisotropy rapidly suppresses the 
"lower" 2:1:1 canted state so that the m=l/3 plateau be- 
comes the third phase under field, above the supersolid. 
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While in the purely nearest-neighbour case this plateau 
emerges continuously connected to canted phases^, in 
this situation it is completely unrelated by symmetry 
to the surrounding phases. Thus this m=l/3 can be 
thought of as an "accident" within the natural progres- 
sion of four-sublattice phases. 

The low-field results at D—0.5 have been analysed 
elsewhere^ and essentially reproduce the Z?=0.25 case. 
At T=0 the first-order phase transition from the super- 
solid into the m=l/3 plateau occurs for a mean-field field 
of 

h =(D + 8( Ji + J 2 ))/3 



- 4/3V2(8J 2 -D- Ji)(D - 2(Ji + J 3 )). (24) 

To study this plateau we use the S z — S z structure fac- 
tor associated with the three-sublattice ordering vectors 
{q K }, cf. Fig. [4) 



^(qK) = (E|^E 5 ? e " iqK ' r 



(25) 



At finite temperature the m =1/3 plateau is con- 
nected through first-order transitions to the m =1/2 



plateau [Fig. 25 , supersolid and 2:1:1 state phases. These 



low-temperature, strongly first-order transitions between 
phases with different symmetries are very difficult to sim- 
ulate, even employing the combined parallel tempering 
and over-relaxation procedure. 

For this value of anisotropy the 2:1:1 canted state sur- 
vives above the m=l/3 plateau for a T=0 field range of 
3.25</i<3.60. While at high field the tendency of the 
m =1/2 plateau to cant through the closing of a spin- 
gap wins, at lower field the m=l/3 plateau is favoured 
entropically against the canted state. Because these two 
transitions lie very close to each other at low tempera- 
ture, it is much more difficult to find them with the same 
accuracy as used in the rest of the phase diagram. 

The nature of the transitions into the supersolid phases 
is unaffected by the increase of anisotropy. However the 



magnetisation jumps [Fig. 24 'b)] have become more dra- 
matic. Hence the points with 8m z = connected to the 
m = 1/2 plateau, at h ~ 1.4 and h 3.6, become more 
visible. 

The endpoint of the inner plateau-plateau transition is 
very close to the critical point of the supersolid transition 
at h « 2,T « 0.3, in a shape curiously reminiscent of 
Fig. [21] The hysteresis associated with these strongly 
first-order transitions, especially the plateau-plateau one, 
makes a more in-depth analysis of this region difficult to 
perform. Nevertheless, as we shall see below, this feature 
is not robust for other values of anisotropy. 

Another possibly interesting feature is that at an 
higher field of ftps 1.90 the plateau-plateau transition also 
supports a point for which the magnetisation discontinu- 

As in the previous case, 



ity apparently vanishes [Fig. 25 




we find no evidence of the emergence of higher symme- 
tries or other exotic behaviour at this point. 



FIG. 25: (color online) Double phase transition from the para- 
magnet to m=l/2-plateau and then into the m=l/3-plateau at 
D = 0.5, h = 1.89. Temperature dependence of the structure fac- 
tor measured at momenta corresponding to four-sublattice or- 
dering iS zz (qM) (a), and (b) momenta corresponding to three- 
sublattice ordering S zz (q«) (b). Both paramagnet-m=l/2-platcau 
and plateau-plateau phase transitions are first-order, as seen in the 
energy histograms in insets to (a) and (b), however the magnetisa- 
tion (c) does not exhibits the characteristic jump across the inner 
transition. 



VII. INTERMEDIATE TO STRONG 
ANISOTROPY, £>=0.675 

By further increasing the anisotropy strength we wish 
to verify the stability of the D = 0.5 conclusions, espe- 
cially in the region surrounding the m=l/3 plateau. We 
expect the collinear phases to gain ground and eventu- 
ally suppress the supersolid states. Therefore we choose 
to study D — 0.675 as a representative value. We note 
that the set of parameters with this value of D gives the 
bet fit to the magnetic excitations in AgNiO^pl. 

The general shape of the phase diagram [Fig. 26 
is similar to the previous ones. The transition fields 
have become less dependent on temperature and the 
phases themselves are more widely separated in the 
magnetisation-temperature phase diagram, a reflection 
of the stronger anisotropy. 

The main difference when compared to lower 
anisotropy is the complete suppression of the 2:1:1 canted 
phase. This phase ceases to exist at T=0 for _D>0.625 
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FIG. 26: (color online) Magnetic phases of a layered triangular-lattice antiferromagnet with Ji=l, J2=0.15, J±=— 0.15 and 
easy-axis anisotropy D=0.675. (a) Phase diagram as a function of temperature and magnetic field. Open symbols on the h-axis show 
transitions obtained in mean-field theory. Phase boundaries at finite temperature are obtained from Monte Carlo simulation for a cluster 
of 24x24x8 spins, and determined by peaks in the relevant order parameter susceptibilities. All phase transitions are first-order, except 
where shown with a dashed line. Thick purple dashed line is obtained through a Landau expansion for the supersolid transition, (b) Phase 
diagram as a function of temperature and magnetization. Solid lines running left-right show cuts at constant magnetic field h taken from 
simulations. The coexistence regions associated with first order phase transitions are coloured white. Thick purple dashed lines show 
phase boundaries obtained through a low-T expansion. Further increasing anisotropy leads to suppression of supersolid phases. 



resulting in a first-order plateau-plateau transition for 
h = 6(Ji — 3J2) = 3.3, a value which holds up to the 
Ising limilP^. 

The points where both the collinear-supersolid and the 
plateau-plateau transitions terminate on the Z3 transi- 
tion line are now well separated. They occur approxi- 
mately at the same field /iwl.85 but at clearly different 
temperatures, T«0.19 and T«0.34 respectively. Hence 
this implies a direct transition between the stripe phase 
and the m=l/3 plateau as a function of field. We thus in- 
terpret the aforementioned closeness of these two points 
at D—0.5 as purely accidental for that specific value 
anisotropy, and not a robust feature of the model. 

The previously identified points where the magnetisa- 
tion jumps disappear at the plateau-stripe critical point 
and at both plateau-plateau and paramagnet-plateau 
transitions are still clearly visible. We have checked that 
these features hold for several values of intermediate tD 
and thus can be said to be characteristic of the m=l/2 
plateau phase. 



all temperatures in the magnetisation-temperature phase 
diagram. The increased anisotropy favours the collincar 
state, driving the transition first order, instead of the 
symmetry- allowed second-order transition observed be- 
fore. At T=0 the mean-field transition field is h=6.57, 
which is lower than the field predicted by the closing of 
the spin-wave gap Eq. ([22]), /i c (D=0.675)=6.72. This 
crossover to a first-order transition happens at £>«0.525. 



The simulation results agree with this interpretation, 
showing for all temperatures the hallmarks of a first order 
transition, jumps in order parameters and magnetisation 
and double-peaked energy distribution [Fig. 



27 . The 



nature of the paramagnet-canted transition is unchanged, 
displaying a critical endpoint at h~7.2 and T«0.06 where 
fluctuations at higher temperatures drive the transition 
first-order. 



A relevant difference is that the 3:1 canted phase has 
now become separate from the parent m—1/2 plateau for 
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FIG. 27: (color online) Double first-order transition at 
D = 0.675, h = 6.4 from the paramagnet into 3:1 canted state, and 
from 3:1 canted state into the m =1/2 plateau. Temperature de- 
pendence of the structure factor 5 zz (qivi) (a), U(l) order param- 
eter (b) and magnetic susceptibility (c). The first-order charac- 
ter of both transitions is clear in the order parameter jumps and 
double-peaked energy distributions close to the critical tempera- 
tures in insets to (a) and (b). The lower transition is continuous 
for smaller D, e.g. D=0.25 [Fig.jTs], 



VIII. STRONG ANISOTROPY, 73=1.5 

Increasing the anisotropy strength eventually leads to 
the suppression of all non-collinear phases. At T=0, the 
supersolid state is squeezed out of existence at 7?>0.9 
for J 2 =0.15 by the collinear phases which surround it. 
For anisotropy values of 73>>1.15 the 3:1 canted phase 
disappears and there is only a direct plateau-saturation 
transition, i.e. a magnetisation jump from m=l/2 to 
m=l with increasing field. We thus focus on 75=1.5 as a 
representative value. 

The ordered observed phases are those collinear phases 
which can be accommodated in a three- or a four-site 
unit cell, i.e. the collinear stripe phase at low field and 
then the m=l/3- and m=l/2 magnetisation plateaux. 
At T=0 the transition fields are independent of the 
anisotropy value and agree with the established results 
for the Ising model. All phase transitions are first order 
for all temperatures and fields. 

Nevertheless, even in this strong-anisotropy regime the 
features above identified previously as robust character- 



istics of the model are still clearly visible. These include 
the persistence of the first-order transition between the 
stripe-collinear and m=l/2 phases, and the points with 
Am=0 in both the paramagnet-plateau and plateau- 
plateau transitions. 



IX. THE ISING LIMIT, D=oo 

It is instructive to compare the strong anisotropy re- 
sults with the pure Ising limit, D — ¥ oo, previously inves- 
tigated at T=0 [Ref. US] . The most obvious difference is 
that in the Ising limit the m=l/3 plateau is now directly 
connected to the paramagnetic state. The two values 
of field where the magnetisation jump at the first-order 
paramagnet-m=l/2 plateau transition are still present, 
but one of them now arises where both plateaux and the 
paramagnet meet. The inner transition fields all show 
very little dependence with temperature, indicating the 
validity of the mean-field picture and that all three phases 
possess roughly the same amount of entropy. As a final 
comment, we observe that the zero field Neel tempera- 
ture is a monotonically function of anisotropy. Empiri- 
cally, this can be written as T^{D) = Tjv(oo)— 2 2i 8 (4)+D ' 
until saturating at T^(oo) — 1.127(3) in the Ising limit. 

X. DISCUSSION OF RESULTS 



TABLE I: Catalogue of phases and associated phase diagrams 
for different values of D. 
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In this paper, we have explored the magnetic-field 
properties of a frustrated easy-axis Heisenberg model, 
originally introduced to explain the magnetic properties 
of the hexagonal delafossite 2H-AgNi02 [42] . We have 
used extensive Monte Carlo simulations, combined with 
Landau theory and spin-wave analysis, to obtain a set 
of magnetic phase diagrams for a set of values of single- 
ion anisotropy, ranging from the Heisenberg [73=0] to the 
Ising [D=oo] limits. The model exhibits a spectacularly 
rich set of phase diagrams, with a wide variety of com- 
peting collinear and non-collinear states, some of which 
are magnetic supersolids in the sense of Matsuda and 
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FIG. 28: (color online) (a) Magnetic phases of a layered triangular-lattice antiferromagnet with Ji=l, J2=0.15, J±=— 0.15 and 
strong easy-axis anisotropy D=1.5. (a) Phase diagram as a function of temperature and magnetic field. Open symbols on the h-axis 
show transitions obtained in mean-field theory. Phase boundaries at finite temperature are obtained from Monte Carlo simulation for a 
cluster of 24x24x8 spins, and determined by peaks in the relevant order parameter susceptibilities. All phase transitions are first-order, 
(b) Phase diagram as a function of temperature and magnetization. Solid lines running left-right show cuts at constant magnetic field h 
taken from simulations. The coexistence regions associated with first order phase transitions are coloured white. Thick purple dashed lines 
show magnetisation curves obtained through a low-T expansion. At high anisotropy non-collinear phases have disappeared and Ising-like 
physics is recovered. 



Tsuneto or Liu and Fisher. A catalogue of these phases, 
together with the values of D for which they occur, is 
given in Table |TJ 

The key to understanding the properties of the model 
for small values of anisotropy lies in the Heisenberg limit, 
D = [Fig. 6j. Here the system breaks up into two, 
decoupled, sublattices [Fig[7], and the ordered states seen 
at finite temperature are selected by considerations of 
entropy, rather than energy. However the difference in 
entropy between these phases is very small, and so any 
finite magnetic anisotropy can stabilize new phases at low 
temperature, as observed for D = 0.02 [Fig. [To]. These 
new phases include the collinear stripe state observed in 
AgNiC>2 [39], and the novel magnetic supersolid [Fig[2|c)] 
introduced in [41] . 

That such a supersolid should be stabilized by easy-axis 
anisotropy might at first seem surprising, since easy-axis 
anisotropy naturally suppresses the in-plane magnetiza- 
tion which lends a magnetic supersolid its "superfluid" 
character. In this case, however, the supersolid is driven 
by the balance of competing exchange interactions and 
anisotropy, and not by any delicate order-from-disorder 
effect. This energetic origin of the supersolid makes it 
very robust — it is present at low temperatures for a 



wide range of values of anisotropy < D < 0.9, as il- 
lustrated in the phase diagrams for D=0.02 [Fig. 
£>=0.25 [Fig.p], D=0.5 [Fig. [24] and L>=0.675 [Fig. 



The transition from the collinear "stripe" state 
[Fig. |2ja)] into the magnetic supersolid, as a function of 
magnetic field, can be viewed as the condensation of spin 
wave excitations with finite momentum. This sets the 
new state apart from other, known examples of magnetic 
supersolids, where the mode which condenses has zero 
momentum. As a consequence the supersolid does not 
simply interpolate between collinear stripe (solid) and 
spin-flop (superfluid) states, or terminate in a tctracriti- 
cal point, as envisaged by Liu and Fisher^. 

Instead, at higher temperatures, the supersolid is 
"squeezed out" of the magnetic phase diagram by com- 
peting phases with higher entropy - - the collinear 
"stripe" phase from which it descends, its canted "spin- 
flop" analogue [Z?=0.02; Fig. [To] and, at higher values 
of D, collinear m = 1/3 [L»=0.25; Fig. p] and L>=0.5; 



Fig. 24 and m = 1/2 magnetization plateaux [£>=0.675; 
Fig. 26 . We anticipate that increasing the second- 



neighbour interaction J-2 (fixed in these simulations at 
J 2 = 0.15Ji), will make the magnetic supersolid more 
robust at higher temperatures. 
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FIG. 29: (color online) (a) Magnetic phases of a layered triangular-lattice Ising antiferromagnet with J±=l, J2=0.15, J±=— 0.15. (a) 
Phase diagram as a function of temperature and magnetic field. Open symbols on the h-axis show transitions obtained in mean-field theory. 
Phase boundaries at finite temperature are obtained from Monte Carlo simulation for a cluster of 24x24x8 spins, and determined by 
peaks in the relevant order parameter susceptibilities. All phase transitions are first-order, (b) Phase diagram as a function of temperature 
and magnetization. Solid lines running left-right show cuts at constant magnetic field h taken from simulations. The coexistence regions 
associated with first order phase transitions are coloured white. 



A similar story is repeated with the two other mag- 
netic supersolid phases found at higher values of magnetic 
fields : an unusual 2:1:1 canted state [Fig. [2]^e)], and a 
3:1 canted state [Fig. [2^g)]. The 2:1:1 canted state inter- 
polates between the magnetic supersolid and a collincar 
rn=l/2 plateau. Like the supersolid, it owes its stability 
to a balance of competing interactions, and is displaced 
at high temperatures by the collinear m=l/2 plateau, 
which has superior entropy. It not only exhibits a sub- 
stantial staggered magnetization in the S x -S y plane, but 
also a finite (if small) uniform magnetization m±, per- 
pendicular to the magnetic field. In turn, the 3:1 canted 
phase interpolates between the collinear m=l/2 plateau 
and saturation. Like the 2:1:1 canted phase it exhibits 
a small magnetization m± in the S x -S y plane. It is the 
only supersolid phase present for D = (where m± = 0) , 
[Fig. [6] , and the only supersolid phase favoured by con- 
siderations of entropy alone. For this reason it is also the 
only supersolid connected to the paramagnetic phase at 
high temperatures. 

The phase transitions which link these different phases 
also evolve as a function of D. Each of the canted phases 
is connected to neighbouring collinear phases by a soft 
spin-wave mode within the collinear state, and so these 
transitions can be continuous. However phase transitions 
at high temperatures are generically first order, and be- 



come more strongly so as anisotropy is increased, as ob- 
served by e.g. comparing the phase diagrams for D = 0.5 
24], D = 0.675 [Fig. [26] and D = 1.5 [Fig.p] 



[Fig. 



This 

trend is seem most clearly in phase diagrams plotted as 
a function of temperature and magnetization, which be- 
gin to exhibit large regions of phase coexistence as D is 
increased. 

The anaysis in this paper clarifies the origin of the 
novel magnetic supersolid introduced in |41j , explores its 
unusual properties at finite temperature, and confirms 
its robustness against other competing phases for a wide 
range of parameters. However a number of interesting 
questions remain. One of these is, how would the system 
respond for magnetic field not aligned with the easy axis ? 
An obvious limiting case is a field perpendicular to the 
easy axis. In this case, the two sublattices of the collincar 
stripe phase can respond to field simply by canting, and 
(zero temperature) mean-field calculation indicate that a 
canted stripe phase interpolates to saturation for all finite 
D. The second obvious limit is that of a field at a small 
angle to the easy axis. In this case, the states described in 
this paper will generally survive, but the symmetries they 
break will be modified by the presence of a component 
of magnetic field in the S x -S y plane. This in turn will 
lead to a modification of the phase transitions between 
them. The limit of small angle is also of relevance to 
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experiments which measure magnetization through mag- 
netic torque, discussed in the context of AgNiC>2 above. 
Away from either of these limits, for intermediate angle, 
the story is not so simple — the balance of energy and 
entropy can favour many different states. This remains 
an interesting topic for future study. 

Another interesting open question is the study of the 
quantum mechanical effects in this spin model. The suc- 
cession of continuous phase transitions observed in our 
model for moderate values of anisotropy, raises the in- 
triguing possibility that a quantum magnet like AgNiC>2 
could exhibit a series of quantum phase transitions as 
a function of magnetic field. Concentrating on the su- 
persolid phases, two-dimensional quantum models with 
broadly similar Hamiltonians support su persol ids both at 
T=0 [Ref. I33H35] and finite temperature^ 21 . However, 
recent calculations suggest that the supersolid phase 
in the original three-dimensional model considered by 
Liu and Fisher may not survive quantum fluctuations^. 
Each case therefore needs to be considered on its own 
merits. One simple, phenomenological, way to incorpo- 
rate quantum fluctuations is through the addition an ef- 
fective biquadratic interaction term, mimicking the se- 
lection of collinear states by quantum fluctuations^!. We 
have checked within classical Monte Carlo simulations 
that such a term does not change the nature of the tran- 
sition into the supersolid phase. On the strength of this, 
and on general grounds, we anticipate that the supersolid 
introduced in |41) will prove robust against quantum ef- 
fects. However this remains to be tested. 



XI. APPLICATION TO EXPERIMENT 

The advent of high field facilities offering static fields 
of up to 45T, and pulsed fields of up to 600T, has made 
it possible to explore the high-field properties of a wide 
range of magnets for the first time. Frustrated mag- 
nets are foremost among these, with model systems like 
TICUCI3 providing a perfect opportunity to study how 
the Bose-Einstein condensation of magnons gives rise to 
a magnetic superfluid [Ref. [3] . Generically, these systems 
exhibit some degree of magnetic anisotropy, either at the 
level of a single ion, or in their exchange interactions. 

Triangular lattice antiferromagnets are no exception. 
Examples with weak easy-axis anisotropy (relative to ex- 
change interactions) inc lude the quasi- two dimensional 
halides VBr2 and VClsF^H which contain antiferro- 
magnetically coupled spin-3/2 V 2+ ions on a trian- 
gular lattice. The insulating oxide Rb4Mn(Mo04)ap21 
and the multiferroic material RbFe(Mo04)2^"^, are 
well-described by a spin-5/2 nearest-neighbour Heisen- 
berg antiferromagnet on a triangular lattice with mod- 
erate easy-axis anisotropy (Rb4Mn(Mo04)3), or easy- 
plane anisotropy (RbFe(Mo04)2). The multiferroic ma- 
terial KFe(Mo04)2 has similar underlying chemistry to 
RbFe(Mo0 4 ) 2 , but a distortion of the triangular lattice 
leads to a somewhat more complex phenomenology^. 



Much stronger single-ion anisotropy, of order of the ex- 
change interactio ns, is found in the hexagonal Ni halides 
CsNiCl 3 , CsNiBr^and in the Cr oxide LiCrOj^ 11 , where 
triangularly coordinated 5=3/2 spins order on the mag- 
netic easy-axis. The insulating oxide 

Further examples are found in the delafossite fam- 
ily, including 5=5/2 CuFeO^, 5=3/2 PdCrO^ and 
5=1 AgNiOj^EH with easy-axis anisotropy, and 5=3/2 
CuCr02- ? with easy-plane anisotropy. Notably, all of 
these delafossite materials require a frustrated model 
such as Eq. ([I]), with second (or further) neighbour in- 
teractions, to accurately describe their magnetic proper- 
ties. Here we have concentrate on AgNi02, where pow- 
der neutron scattering studies provide clear evidence for 
a collinear, stripe ground state in zero magnetic fiekP^, 
as would be expected for the frustrated, easy-axis model 
considered in this paper with moderate J2 and D. 

At present, single crystals of AgNi02 are too small 
for inelastic neutron scattering experiments to be per- 
formed on them. However angle-integrated spin-wave 
spectra can be measured in powder samples. These 
show clear evidence of a spin-gap, and spectra are well- 
described by the present model [Eq. Q], with parame- 
ters Ji=1.32meV, J 2 =0.20meV, L>=1.78meV, Jj_=-0.14 
meV0H. For spin 5=1, the effective value of D seen in 
quantum spin wave spectra is renormalized by a factor 
two relative to its classical valued and so comparison 
should be made with classical Monte Carlo simulation 
results for D—0.675xJ 1 . Thus, under the (strong) as- 
sumption that this simple spin model provides an ade- 
quate description of AgNi02 in high magnetic field, the 
field-temperature phase diagram for AgNi02 should be 
of the form shown in Fig. 26 [£>=0.675]. Given the 
simplifications inherent in using any spin model to de- 
scribe a metal, the critical fields and transitions temper- 
atures predicted here should be approached with some 
caution, are unlikely to be quantitatively correct. How- 
ever our Monte Carlo simulation results should provide 
a reasonable first guide to the different phases occuring 
in AgNi02 in magnetic field, and the nature of the phase 
transitions between them. 

The actual phase transitions which occur in AgNi02 
for fields of up to 40T, have been studied through meau- 
rements of magnetic torque f = m x h, heat capacity 
and electrical transport"™ These experiments indicate 
that a continuous (or very weakly first order) phase tran- 
sition out of the collinear stripe state occurs at about 
12.5T (at a base temperature of 1.5K). This phase tran- 
sition is accompanied by a weak anomaly in the specific 
heat. And, crucially, it occurs at a critical field which 
increases steadily with temperature, indicating that the 
collinear stripe phase has a higher entropy than its suc- 
cessor. These are exactly the characteristics of the phase 
transition from the collinear stripe phase into the super- 
solid for the frustrated easy-axis model Eq. ([I]). 

In Ref. [41] we therefore proposed that a novel mag- 
netic supersolid is realised in AgNi02 for magnetic fields 
greater than 12. 5T. We further suggested that torque 
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was a good tool for distinguishing phase transitions in 
easy-axis magnets, exhibiting changes in sign as well as in 
magnitude, and made explicit predictions for the torque 
signature of the proposed magnetic supersolicPS. It is 
also interesting to note that torque clearly distinguishes 
those phases, such as the 2:1:1 canted phase and the 3:1 
canted phase above, which have a finite magnetization 
m± in the S x -S y plane. 

The results of this paper confirm that the scenario pre- 
sented in [4Tj is robust across a wide range of parame- 
ter space, and does not require any special assumptions 
or fine-tuning of the model. In this context it would 
be interesting to look for evidence of the spin gap clos- 
ing at the putative supersolid transition in AgNiC>2, ei- 
ther from powder neutron scattering, or NMR relaxation 
rates. However it is important also to remember this 
model is a gross simplification of the physics of AgNi02, 
which contains itinerant electrons as well as local mo- 
ments. The development of a more realistic model for 
AgNiC>2, which takes these itinerant electrons into ac- 
count, remains an important avenue for future study. 



XII. CONCLUSIONS 

In this paper we have studied a realistic, three- 
dimensional spin model motivated by the hexagonal de- 
lafossite 2H-AgNi02. We have obtained the full magnetic 
phase diagram of this model as a function of temperature 
and magnetic field for values of easy-axis anisotropy D 
ranging from the Heisenberg (D=0) to the Ising (D=oo) 
limits. We uncovered a rich variety of different magnetic 
phases, including several phases which are magnetic su- 
persolids (in the sense of Matsuda and Tstuneto or Liu 
and Fisher), one of which may already have been ob- 
served in AgNiOsp^- We explored how this particular 
supersolid, first introduced in [41], arises through the 
closing of a gap in the spin-wave spectrum, and how 
it competes with neighbouring phases as the easy-axis 
anisotropy is increased. This novel phase was shown 
to have qualitatively different finite-temperature prop- 
erties from any previously studied magnetic supersolid, 
and to be remarkably robust against changes in parame- 
ters. These results suggest that magnetic supersolids in 
frustrated systems can have a richer phonomenolgy, and 
be far more robust, than previously supposed. 
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Appendix A: Calculation of spin stiffness 

Spin stiffness is defined as a generalised elasticity co- 
efficient representing the free-energy cost of applying a 
twist to the boundary conditions with gradient S/fi along a 
given direction e. Since there is a direct mapping between 
spin stiffness and superfluid density^, a phase with non- 
vanishing ps is said to have superfluid character. In finite 
magnetic field it is sufficient to consider twists around S z , 
i.e in the S x -S v plane. For a single pair of spins we write 

Si • = COS (>; - <j>j + 5<j>.e.(Ti - Tj)), (Al) 

where each spin is expressed in polar coordinates 

Si = (cos(f>ism9i,sm(j)ism9i,cos0i). (A2) 

Spin stiffness is then given by the second derivative of 
the free energy with respect to the twisting angle 5<p 



_ d 2 T _ i d 2 U \ 1 // &H V 

Ps[Gl ~ WW s*=o~ \ WW' H=o T\\d(8$)) 



(A3) 



Applying Eq. ( A3 ) to both first- and second-neighbour 



interactions and normalising per unit area (i.e. per spin), 
spin stiffness readfp2E21 



Ps(e) 



V3NT 



^| v /j 1 ^(e.r y ) 2 .S^S^+J 2 ^(e.r y ) 2 .S^S 
/{j 1 ^(e.r iJ ).S^xS|+J 2 ^(e.r y ).S^xS|} 2 



(A4) 



where ry=rj— Yj. Since the parallel tempering method 
restores the full lattice symmetries, we average ps 
over the three symmetric directions in the lattice 
e=(e x , e y )= {(1, 0), (1/2, V3/2), (-1/2, y/3/2)}. 



Appendix B: Low temperature expansion and 
Landau theory for supersolid transition 

The low-temperature excitations of a classical spin 
model of the form Eq. ([I]) are frozen spin waves, in con- 
trast with the dynamic spin waves found in quantum 
magnets. These excitations have a dispersion which can 
be calculated through an expansion of each spin fluctu- 
ations around its preferred T—0 configuration. In the 
sufficiently generic case of two spins in a common plane 
canted by an angle 29, the Heisenberg interaction is ap- 
proximated by 



Si-Sj « cos 26»- -(xj 
+ XiX4 cos 26 - 



ViVj, 



y 2 ) cos 29 



(Bl) 
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where x% and yi denote small fluctuations about the or- 
dered state in the local Sf and Sf directions respectively. 
The rotation into a common frame is arbitrarily chosen 
to be performed around the S v direction. 

The resulting low-temperature Hamiltonian is solved 
by Fourier transform, with two spin- wave modes e a (k) 
per n sublattice, i.e. 



tL — rlo 



^ 2n 



N/n 



2 ^ ^ 

a keM B 2 



( (k)x c 



(B2) 



We now illustrate this procedure for the supersolid 
state with four different sublattices : while two are 
collinear with S z in the positive direction the other two 
are canted by an angle of 8 from S z in opposite directions. 
After writing down the interactions between sublattices 
in matrix form, the spin-wave dispersion modes can be 
obtained by diagonalisation of 



M£ = 2 
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where eg— cos 9 and the other coefficients are given by 
He = (Ji + J 2 )(2c e - 1) + J±( 7z (k) - 1) + D ■ 
W k = (Ji + J 2 )(2c e - c 2B ) + J±(7,(k) - 1) - '-co, 



h 
2' 



X k = Ji 7 r(k) + J 272 ^(k), 
y k = -Ji7i + (k)- J 2 7 2 ~(k), 
Z k = -J l7l -(k)- J 272 + (k). 



(B5) 



The lattice structure factors are 
7 5 4B (k)=cosfc a; , 72 4B (k)=cos v^fcy, 7z (k)=cosfc z , 
7 j t (k)=cos((fc 2; ± v / 3fc 2/ )/2) and 7 ^ (k)=cos((3fc x ± 
\/3ky)/2). The four-sublattice structure results in 8 
dispersing modes e(k). The thermodynamics at low-T 
are accessed via the free energy, 



T 

N 



Eo 
N 



2n N/n 

- TlnT+-^Yl ^lne Q (k)+0(T 2 ). (B6) 
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The free energy of the supersolid phase near its phase 
transition can be obtained as a Landau expansion in pow- 
ers of the order parameter measuring the U(l) broken 
symmetry 



(7(1) 



(B7) 



The U(l) order parameter [Eq. ([5])] is by construction 
proportional to the canting angle 0, as represented in the 
cartoon Fig. [2jc) , and therefore near the transition 



O u{1) = 2sin0 ~ 0. 



(B8) 



The coefficient a collects the quadratic contributions 
from the spin wave modes 



a = hssoiT) 
hsso (T) = 2D + 2T^{^J2 ln e « ( k ) } 



h, 

-)2 



a,k 



0=0 



(B9) 
(BIO) 



where the critical field is calculated near T=0. The 
prefactor for the temperature term is obtained by nu- 
merically integrating the derivatives of the eight spin- 
wave branches over the four-sublattice magnetic Brillouin 
Zone, corresponding to the finite lattice used in simula- 
tion. The resulting critical fields are plotted as a yellow 
line on the (a) panel of Figs. |10|13[ [24] and [26], yielding a 
very good agreement with simulation, even at relatively 
high temperatures. 
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